Reverse-selective diffusion in nanocomposite membranes 
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The permeability of certain polymer membranes with impenetrable nanoinclusions increases with 
the particle volume fraction (Merkel, et al., Science, 296, 2002). This intriguing observation con- 
tradicts even qualitative expectations based on Maxwell's classical theory of conduction/diffusion 
in composites with homogeneous phases. This letter presents a simple theoretical interpretation 
based on classical models of diffusion and polymer physics. An essential feature of the theory is a 
polymer-segment depletion layer at the inclusion-polymer interface. The accompanying increase in 
free volume leads to a significant increase in the local penetrant diffusivity, which, in turn, increases 
the bulk permeability while exhibiting reverse selectivity. This model captures the observed depen- 
dence of the bulk permeability on the inclusion size and volume fraction, providing a straightforward 
connection between membrane microstructure and performance. 



Polymer membranes facilitate a variety of molecular 
separations. Because their micro-structural scale is com- 
parable to the size of solute molecules, the polymer ar- 
chitecture can be tailored to separate specific gas mix- 
tures. Merkel et al. @ recently showed that incorpo- 
rating nanometer sized inorganic particulates into cer- 
tain amorphous polymer melts increased the membrane 
permeability and selectivity Because of its significant 
technological applications, this discovery has stimulated 
many experimental investigations [3. IsL flTL IT^ . A theory 
that quantifies how the inclusion size and concentration 
affect the permeability and selectivity does not exist. In- 
deed, Merkel et al. highlighted that classical Maxwell-like 
theories fail to describe even the qualitative trends. 

This letter presents a theoretical interpretation based 
on classical models of diffusive transport and polymer 
physics. Despite its simplicity, the model accurately de- 
scribes experimental trends. The principal assumption 
is that there exists a repulsive interaction between the 
nanoinclusions and polymer segments. This hypothesis 
is support, in part, by transmission electron microscopy 
(TEM) images of silica-based nanocomposites [e.<7.,0,|l|, 
which show clear evidence of particle aggregation. 

For steady gas permeation across a membrane with 
thickness L, the diffusive flux is 

\{j)\ = -D e An/L, (1) 

where D e is the effective diffusivity and An is the differ- 
ential concentration of the diffusing solute. Because the 
solute concentrations at the gas-membrane interfaces are 
proportional to the partial pressure p, with a constant of 
proportionality H (Henry's constant), it is customary to 
express the flux in terms of the differential partial pres- 
sure Ap: 

= -D e HAp. (2) 
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The permeability is defined as |(j)|/|Ap| = D e H. In 
this work, the gas solubility is assumed to be indepen- 
dent of solids concentration, so the inclusions influence 
the permeability only through their role in modifying the 
effective diffusivity D e . 

If voids in the polymer are much larger than the solute 
molecules, then the diffusivity in these voids D will be 
much larger than in the polymer D°° . When D/D°° — > 
oo, Maxwell's theory yields an effective diffusivity D e — 
D°°(l + 3<f) v ) + 0{4>1), where <p v is the (small) volume 
fraction of (spherical) voids. Under these conditions, the 
relative enhancement of all effective diffusivities are the 
same, so membrane selectivity is unaltered. 

Merkel et aZ.'s experiments exhibit reverse selectiv- 
ity, meaning that the permeability of larger molecules 
is enhanced more than smaller ones. This necessitates 
molecular-scale perturbations to the polymer microstruc- 
ture, precluding the notion that the permeability is en- 
hanced by a shell of void space between the (impene- 
trable) inclusions and surrounding polymer. Rather, it 
suggests a continuous increase in the free volume when 
approaching the solid-polymer interface. Figurcndepicts 
a single impenetrable sphere embedded in a polymer melt 
where polymer segments are repelled from the solid sur- 
face. This is the microstructural view that underlies the 
theory presented below. 

Cohen and Turnbull's statistical mechanical theory Q 
yields a solute diffusion coefficient 

D = Acxp (-yvm/vf) (3) 

where A and 7 are constants, v m is the minimum free 
volume required for a solute molecule to escape its cage 
of neighboring atoms and, hence, undergo diffusive mi- 
gration, and Vf is the available free volume per volume 
occupying element. In this work, the volume occupying 
elements are the atoms of the polymer melt. For sim- 
plicity, each atom is assumed to occupy, on average, a 
1/3 

volume vq, where v is of the order of a covalent bond 
length (« 1.5 A). By considering the total volume, which 
comprises the sum of free and occupied volume, it follows 
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FIG. 1: An impenetrable nanosphere with radius a ~ 5 nm 
embedded in a polymer matrix with correlation length £ ~ 
1 nm. A repulsive interaction between the inclusion and the 
polymer produces a surface layer (with thickness ~ £) where 
the polymer segment density approaches zero. The local in- 
crease in free volume w/ (within the dash-dotted circle) in- 
creases the penetrant diffusivity, leading to a significant in- 
crease in the bulk permeability when £ > a (see figure 



that 



v f = vo{mi/{v nxp) - 1], 



(4) 



where Co is the average atomic concentration, wi\ is the 
mass of a monomer (repeat unit), n\ is the number of 
atoms per monomer, and p is the (mass) density. 

A tractable analytical expression for the radial polymer 
segment density is obtained from a self-consistent mean- 
field model, with the so-called ground-state approxima- 
tion and a flat interface 0- The influence of curvature 
(finite inclusion radius) is of secondary importance, and 
is therefore neglected here. With a repulsive interaction 
between the polymer segments and the solid, the segment 
concentration is .3] 



3 (r) = Coo tanh 2 [(r - a)/£] + O^/a) 



(5) 



where Coo is the bulk concentration, and £ = 1/ ^/ivcoc is 
the polymer correlation length, with / the segment length 
and v the excluded volume (per segment). Conservation 
of chain contour length and mass requires 



<j> = cl 3 = (p/ mi )hl 2 



(0) 



where l\ is the length of a monomer. Combining 
Eqns. ©-© gives 



D = D°° exp [- 



v m<t>* (<t> ~ <l>oo) I 
(l-r0)(l-0*0oo) J 



(7) 



where <f>* = niVo/(hl 2 ) and = jv m /vo are dimen- 
sionless parameters that reflect the prevailing atomic and 
molecular geometry 

Let us simplify matters by assuming that the bulk 
polymer concentration tpoo = CqqI 3 ~ 1 and « « I 3 , so 



I 2 w mi/(lip), £ ~ l/y/3, and 4>* ~ niv p/m\. Under 
these conditions, 1 — 0* is the fractional free volume of 
the bulk polymer, and the (maximum) diffusivity at the 
solid-polymer interface becomes 



D(r = a)«Z>~expK07(l-0*)]. 



(8) 



Clearly, the diffusivity at the interface can be much larger 
than in the bulk when (j>* < 1 and v* n < 10. Note that 
reverse selectivity prevails, because the diffusivity de- 
creases continuously with increasing radial distance from 
the solid-polymer interface. 

Consider, for example, bulk poly(p-trimethylsilyl 
styrene) (PTMSS), for which m = 28, m 1 = 180.3 g/mol, 
h w 3 A (two C-C bonds) and p w 965 kg/m 3 0. With 
</>«, = 1, Z « 1.02 nm, f w 0.59 nm and 0* = 0.090w 
(with iJg^ 3 measured in A). To determine an appropriate 
value for Vo, which, in principle, should not vary sig- 
nificantly from one polymer to another, let us adopt a 
reported value of the fractional free volume 1 — <f>* = 



0.191 |2| . This provides w 2.0 A and v 1 / 3 w 1.3 A, 
which are both of the expected magnitude. 

With a statistically homogeneous microstructure, the 
average diffusive flux can be expressed as a volume aver- 
age 



(j) = v- 



jdV, 



(9) 



where the integration extends over the discrete and con- 
tinuous phases of a representative elementary volume V. 
The local diffusive flux is j = —D'Vn (Fick's first law) 
where n is the solute concentration. Under steady con- 
ditions, conservation demands 



V • (DVn) = 0, 



(10) 



where D(r) is given by Eqns. (JSJ and Q. When the 
volume fraction <j> p = n p (4/3)7ra 3 <C 1, the average flux 
is H 

(j) «-^ oo (Vn)+30 p (S/a 3 )^ oo (Vn) + O(02), (n) 

where (Vn) = HAp/L is the average concentration gra- 
dient and B is the dipole strength, i.e., 



(Vn) ■ r + BCVn) ■ rr 



oo. 



(12) 



In this work, B is obtained from an efficient numerical 
solution of Eqn. (|10JI that satisfies Eqn. (I12|) as r — > oo 
with a no-flux boundary condition at r = a|l3|. 

The results are calculated with (p^ = 1 and <fi* = 0.8. 
The two remaining independent (dimcnsionlcss) param- 
eters are the scaled correlation length £/a and the scaled 
tracer size v* — ^yvm/vo. Here, the principal influence 
of v* m is to set the diffusion coefficient at the interface: 
D(r = a)/D°°=exp(4v* m ). 

As expected, all values of v* n > increase the effective 
diffusivity, because the available free volume increases 
with respect to the unperturbed (bulk) polymer. Also, 
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FIG. 2: The scaled diffusivity increment A = (D e /D x - 
1)/4> P versus the scaled polymer correlation length £/a with 
— jv m /vo = 0.1, 0.2, 0.4 and 1.0 (increasing upward), and 
cj>* — 0.8 (solid lines). The dashed lines are exact solutions 
(Eqn. 1141 for impenetrable inclusions with a uniform coating 
of thickness £ where D(a < r < a + = exp (4v* n ). 



FIG. 3: The scaled effective diffusivity D e /D°° versus the 
scaled inclusion radius a/£ with v* m = yVm/vo = 0.1, 0.2, 
0.4 and 1.0 (increasing upward), (f>* = 0.8, and <f> p = 0.13: 
0(4> P ) (solid lines); O(cj)p) (dashed lines). The circles are 
experimental measurements of the permeability enhancement 
from || (reported radii scaled with £ = 0.8 nm). 



because the relative increase in diffusivity depends expo- 
nentially on the penetrant size (Eqn. ©), situations with 
multiple diffusing species exhibit reverse selectivity @ . 

Figure presents the most general predictions of the 
theory. The diffusivity increment, defined as 



A = -3B/a 3 = (D e /D°° - l)/<j> p + 0(<j> p ), 



(13) 



is plotted as a function of the scaled correlation length 
£/a for various values of (solid lines). These calcu- 
lations are compared with the exact solution for an im- 
penetrable sphere with a uniform coating (dashed lines): 



A = -3(l + f/a) 3 (1 - /3)[1+a(1+W]+3/3 (14) 

Here, ft = D/D°°, with D the diffusivity inside the coat- 
ing (a < r < a + £), £ is the coating thickness, and 
a = 2 [l4|. The figure demonstrates that the effec- 
tive diffusivity from Eqn. JSJ is much more sensitive to 
the layer characteristics, particularly when the layers are 
thick (£/a > 1). 

To compare the theory with the available experimental 
data, figurcs[3]and0]show the relative effective diffusivity 
D e /D°° and measured values of the relative permeabil- 
ity D e H/(D 00 H 00 ) (circles), with cj) p = 0.13. The 0(d> p ) 
theory (solid lines) neglects particle interactions, so the 
values are as given by Eqn. Hll|) . The dashed lines are 
an 0(<p 2 ) theory, which has elements of a self-consistent 
mean-field approximation with an explicit correction for 
interactions between pairs of particles in a statistically 



homogeneous dispersion [g. Note that the (single parti- 
cle) dipole strength and particle concentration both affect 
two-body interactions. 

Some of the experimental scatter in figure [S] may be 
attributed to the variety of filler particles (all embed- 
ded in poly (4- methyl- 2-pentyne) (PMP)) and, possibly, 
different penetrant species and degrees of particle ag- 
gregation Nevertheless, with the foregoing approx- 
imations, the correlation length inferred by the fit is 
£ = 0.8 nm (with <fi* = 0.8) and, hence, the seg- 
ment length I f=a ~ 1-4 nm. The repeating unit 
of PMP comprises two C-C bonds, so there are about 
4.5 monomer units per statistical segment, which is rea- 
sonable because polymers with much less "bulky" side 
groups (e.g., poly(oxyethylene)) have fewer than two re- 
peat units per statistical segment [e.g. flfl pg. 168]. 

Figure 0] shows how the effective diffusivity (with 
a = 6.5 nm) increases with the solid volume frac- 
tion. The theory is presented with a correlation length 
£ = 0.8 nm, which follows from the data in figure |3 so 
£/a = 0.8/6.5 w 0.123. Despite the experimental data 
extrapolating to a value of D e < D°° as <f) p — > 0, the 
theoretical and experimental trends are in good agree- 
ment. For reference, the dash-dotted line is Maxwell's 
self-consistent theory for impenetrable inclusions in an 
unperturbed polymer matrix. It is remarkable, perhaps, 
that a perturbation extending only £ w 0.8 nm from the 
inclusion surfaces can have such a significant influence on 
the bulk permeability. 

Theoretical predictions of the effective diffusivity have 
been compared with experimental measurements of the 
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relative permeability. Therefore, the solubility was as- 
sumed independent of the inclusion concentration. Fur- 
ther, a mean-field description of the polymer segment 
density was adopted, with a flat, repulsive interface where 
the segment concentration vanishes (at r = a). Also, the 
diffusion coefficient was assumed to follow Cohen and 
Turnbull's formula (Eqn. ©), with the occupied vol- 
ume assumed proportional to the atomic number den- 
sity. Finally, theoretical predictions were based on a sta- 
tistically homogeneous dispersion of inclusions, whereas 
the particles are often found to aggregate. Nevertheless, 
the semi-quantitative agreement of theory and experi- 
ment strongly suggests that reverse selectivity in these 
nanocomposite membranes is of classical origin. 



FIG. 4: The scaled effective diffusivity D e /D oc versus the in- 
clusion volume fraction <j> p with v* m = yv m /vo = 1-0, cj>* — 0.8 
and £/a = 0.8/6.5 ~ 0.123: exact 0(<j> P ) theory (solid lines); 
approximate 0((f>p) theory (dashed lines). The circles are ex- 
perimental measurements of the permeability enhancement 
from || (with a = 6.5 nm). The dash-dotted line is Maxwell's 
self-consistent theory for impenetrable inclusions and unper- 
turbed (homogeneous) polymer. 
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